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In this article we have investigated the effect of weak random disorder in the BCS-BEC crossover 
region. The disorder is included in the mean field formalism through NSR theory of superconducting 
fluctuations. A self consistent numerical solution of the coupled equations involving the superfluid 
gap parameter and density as a function of the disorder strength, albeit unaffected in the BCS phase, 
yields a depleted order parameter in the BEC regime and an interesting nonmonotonic behavior of 
the condensate fraction in the vicinity of the unitary region, and a gradual depletion thereafter, as 
the pairing interaction is continuously tuned across the BCS-BEC crossover. The unitary regime 
thus demonstrates a robust paradigm of superfluidity even when the disorder is introduced. To 
support the above feature and shed light on a lingering controversial issue, we have computed the 
behavior of the sound mode across the crossover that distinctly reveals a suppression of the sound 
velocity. We also find the Landau critical velocity that shows similar nonmonotonicity as that of 
the condensate fraction data, thereby supporting a stable superfluid scenario in the unitary limit. 

PACS numbers: 05.30.Jp, 74.20. Fg, 74.40.-fk, 03.75.Lm 



I. INTRODUCTION 

Atomic gases at very low temperature are a unique 
system where one can observe the continuous evolution 
of a fermionic system (BCS type) to a bosonic system 
(BEC type) by changing the inter-atomic interaction by 
means of Fano-Feshbach resonance IJ. The experimen- 
tal advances to address this transition has introduced the 
possibility for studying the so called BCS-BEC crossover 
more closely 2 , 3 . The physics of the crossover focuses 
on the change of s-wave scattering length (1/a) from at- 
tractive to repulsive (from — oo to -l-oo) by tuning an ex- 
ternal magnetic field. In the intermediate region where 
1/a — > 0: a new region emerges, where dimensionless co- 
herence length becomes in the order of unity, is the focal 
point of crossover physics. 

The situation becomes further interesting if one as- 
sumes existence of random disorder in an otherwise very 
clean system. It is well-known that every wavefunction 
is spatially localized when the disorder is introduced in 
one-dimensional case, which is called as Anderson local- 
ization However, it is immensely difficult to directly 
observe the localization in electronic systems on crystal 
lattices, so that one has to take the indirect route of 
conductivity measurement to observe the effects of local- 
ization. It is an intriguing question how the Anderson 
localization modifies the BCS superconductivity. It has 
been found by Anderson 5 that the order parameter is 
unaffected when the disorder is not too strong to give 
rise to the Anderson localization. Here the time reversal 
pairs, rather than the opposite momentum pairs, then 
form the Cooper pairs. 
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The cold atomic system offers great amount of con- 
trollability and allows one to observe macroscopic wave- 
function. Therefore, it was conceived as a very useful 
candidate to visualize the disorder effects directly. Re- 
cently the ultracold Bose gas {^"^Rb and •^^K) enabled us 
to see the localization directly |6l Pfl . Latest experiment 
are conducted in three dimension for both noninteract- 
ing atomic Fermi gas of '^'^K [8j and bose gas of ®''i?6 [S]. 
These experiments has widened the possibility to study 
the crossover in lights of disorder [lOj experimentally. 

The static disorder in Fermi and Bose systems are not 
new issues. A considerable amount of attention has been 
paid to disorder in superconductors [TTHTB] and in Bose 
gas [HHIS] . Of late the interest at unitarity is also gain- 
ing pace PSMT] , but still it failed to address all the ques- 
tions associated with it. It is exciting to envision the 
three dimensional phase diagram involving temperature, 
interaction and disorder through the evolution from BCS 
to BEC superfluid. At the beginning it can be considered 
that the random potentials are independent of the hyper- 
fine states of the atoms but then it can be extended to 
the correlated disorder problem in the crossover region. 

In this article, we present our investigation on BCS- 
BEC crossover with weak uncorrelated disorder at zero 
temperature. Precisely, we show: (i) monotonic deple- 
tion of order parameter, (ii) nonmonotonic nature of con- 
densate fraction, and (iii) suppression of sound velocity, 
as a function of disorder. Further we add a study on 
Landau critical velocity to attempt a qualitative under- 
standing of the nonmonotonic behavior of the condensate 
fraction. Though order parameter is mainly for academic 
interest but the other three quantities are experimentally 
viable [32U34| . In this study, we observe that super- 
fluidity is more robust in the unitary regime, which is 
also consistent with behavior of the coherence length in 
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the crossover regime. The robust nature of superfluid- 
ity has already been pointed out in the context of vortex 
core structure [35], Josephson current [35] and collective 
modes [37]. We also observe a progressive decay of the 
sound velocity in the BEG side with increasing disorder 
possibily occuring due to enhancement of impurity scat- 
tering [3H|. Here we include a systematic study of the 
physical observables and their response to the random 
disorder. 

The recent experiments on ultracold Bose and Fermi 
gas with disorder was carried out with optical speckle 
[5] IS] and quasi periodic optical lattices [7] . Though they 
pose interesting physics to study, but here we consider 
quenched delta correlated disorder which remained sub- 
ject of interest at zero temperature |26) and at finite tem- 
peratures |31j in previous studies. One can visualize the 
situation of impurity driven random potential by con- 
sidering the presence of few heavy atoms (say ^'^K) in 
a homogeneous bath of light ^Li atoms [3T]. Since the 
number of heavy atoms are very limited, one can term it 
as a quasi-homogeneous system. To include the disorder 
effects in the mean field approach, we follow the NSR 
theory [32] of superconducting fluctuations extended to 
broken symmetry state [ID]. Further we take advantage 
of the techniques developed to study condensate fraction 
of clean Fermi gas with Gaussian fluctuations [U |32] • 

We arrange our study in the following way, in Section 
II we present the basic formalism. Section III is ded- 
icated to discuss our results in three parts. First part 
contains the order parameter, second one for the con- 
densate fraction and the third leads to sound mode and 
Landau critical velocity. Finally, in Section IV we draw 
our conclusions. 



II. FORMALISM 

To make our presentation self-contained, we briefly 
summarize the mathematical formalism presented in 
Ref. [26] . To describe the effect of impurity in Fermi su- 
perfiuid in the crossover from BGS to BEG regime one 
needs to start from the real space Hamiltonian in three 
dimension for a s-wave superfluid, 



H(x)=^ci>t(x) 
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J dx'V(x,x')<i>|(x')<I>l(x)$^(x)<i>^(x'), (1) 

where '^^(x) and $(t(x) represents the creation and an- 
nihilation of fermions with mass m and spin state cr at x 
respectively. V(;(x) signifies the (weak) random potential 
and /i is the chemical potential. We set h — 1, where 
h is the Planck constant. The s-wave fermionic interac- 
tion is defined as V(x, x') = —gS{x — x'). The disorder 
potential is modeled as, Vd(x) = J2i9dS{^ — Xi) where 
gd is fermion-impurity coupling constant and x^ are the 
static positions of the quenched disorder. We assume it 



exhibits white noise correlation, that is, {Vd{—q)Vd{q)) = 
i^PSiu^fi. P is the inverse temperature, Vm is the bosonic 
Matsubara frequency and k = '^iff^' ^^^^ describes the 
strength of the impurity potential with rii being the con- 
centration of the impurities. 

The partition function corresponding to Eq. ([T]) can be 
written in the path integral formulation as 



Z= / 2?[$,a>]exp[-5({$}{$})], 



(2) 



1) Eq.p^ can 



where S — dr J dxl^^dr^a + 'H]. By introducing the 
pairing field A(x, r) and by applying the Grassman iden- 

tity iJV[A, A] exp [- / dx /^^ drAA/g] 
be given as 

Z^J 2?[$, $] J 2?[A, A] exp [-5e//], (3) 

where Seff = 5($, $) + 1/gJdx drAA. Follow- 
ing Hubbard-Stratonovich transformation, Eq.(|3]) can be 
written in terms of inverse Nambu propagator as, 



Ze// = J V[A,A]e-'/s^''''^o'i-^^ 

X y 2? [I>$]e^ ■/■o"''^*^"'*, (4) 

where the inverse Nambu propagator Q^^{x, r) is defined 
as, 



A 



(5) 



After integrating out the fermionic fields from Eq.(|4]) we 
are left with the effective action as. 



Seff = I dx. 



dT 



|A(r)| 
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Tv\n{-f3g-\v)} ,(6) 



where r = (x,t). It is important to mention that the 
main contribution in the partition function comes from a 
small fiuctuation, (5A(x, r) = A(x, r) — A where A is the 
homogeneous BGS pairing field. The Green's function in 
can be written as a sum of Green's function in ab- 
sence of disorder (Qq^ = —drl + (V^/2to -|- fi)az + Aa^) 
and a self energy contribution (S = —Vd<Jz + SAa+ + 
6Aa-) which contains the disorder as well as the small 
fiuctuations of the BGS pairing fields. I denotes the iden- 
tity matrix, and ai are the Pauli matrices and ladder 
matrices(i € {x, y, z, -f , — }). 

By expanding the inverse Nambu propagator upto the 
second order one can write the effective action {Seff) in 
Eq.(|6]) as a sum of bosonic action (Sb) and fermionic 
action (Sp). Also it contains an additional term which 
emerges from the linear order of self energy expansion 
(GqS). It is possible to set zero for the linear order if 
we consider Sp is an extremum of Sef / after performing 
all the fermionic Matsubara frequency sums. The con- 
strained condition leads to the BGS gap equation which 
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after appropriate regularization through the s-wave scat- 
tering length reads, 



m 



1 

2l~k 



(7) 



Eq.([7]) suggests that the BCS gap equation does not have 
any contribution from the disorder potential explicitly. 

Now to construct the density equation with usual pre- 
scription of statistical mechanics; the thermodynamic po- 
tential O should be differentiated with respect to the 
chemical potential /i. O can be written as a sum over 
fermionic (flp) and bosonic (yis) thermodynamic poten- 
tials, which implies, 

d \ d 

n = nF+nB = - — (f^F + ^^s) = --.^{Sp + Sb18) 

The well known BCS density equation can be restored 
from Eq.Q if we consider only np, then it yields the 
familiar X]fc(l ~ b^)- However the presence of disorder 
and fluctuation leads to 0- Hence the final mean 

field density equation will be. 



(9) 



The bosonic thermodynamic potential consists of two 
parts. One comes from the thermal contribution and the 
other is due to disorder. Since we are interested in zero 
temperature we neglect the thermal contribution from 
here on. Henceforth the disorder induced thermodynamic 
potential can be written as. 



K 
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J2 J^^M-'M. 



(10) 



where A/" is a doublet which couples disorder with fluc- 
tuation. After performing the fermionic Matsubara fre- 
quency summation over Af, 



2EkEk+q{Ek + Ek+q) 



(11) 



The inverse fluctuation propagator matrix is a 2 x 2 
symmetric matrix whose elements are given by 



2 2 



Ml2 = y^^UkVkUk+qVk+q 
k 

1 



UkUk+q 

Ek — Ek+q iVm + Ek + Ek+q 

1 



iVm + Ek + Ek+q 



Ek — El 



k+q 



(12) 



and M22{<l) — A/(ii(— g). Here the usual BCS nota- 
tions have been used, namely ^fe = k^/2m — iJ., Ek = 
yef+A^, ul = U^ + ^k/Ek) and vl = 1(1 - ^k/Ek). 




FIG. 1. The order parameters A normalized by Fermi energy 
as function of {kpa)'^ for various disorder strengths. The 
shaded area represents the crossover region 



III. RESULTS 



Order Parameter 



Eqs.Q and ([9| are now ready to be solved self con- 
sistently. Our analysis is valid only for the weak disor- 
der. Considering Ref. [3P we can safely assume that the 
disorder is weak if the dimensionless disorder strength 
?](= KTrt^ /kp) < 5 is satisfied. Fig. [l] demonstrates that 
in the BCS limit (l/kpa— > — oo) every order param- 
eters with different disorder strength follows the mean 
field expression of A/ep = 8e^^ exp[— 7r/(2fci?a)] empha- 
sizing the validity of Anderson theorem ^5j. In the BEC 
limit, instead, one can observe a progressive depletion. 
Assymptotically this depletion remains in the order of 
r]/{kpa) [26l[43]. The other mean field quantity, chemi- 
cal potential, does not change much(not shown). We un- 
derstand that it might be an attribute of the fluctuation 
theory ^35j where the correction in fi at the BCS side is 
usually in ©(A^) which is a very small quantity as A — > 
for {kpa)~^ — 0. In the BEC side the correction comes 
through the effective chemical potential of the composite 
bosons. However the BEC chemical potential is domi- 
nated by the binding energy and it turns out quite large 
compared to the effective chemical potential. The exter- 
nal potential Vdix) induced by the disorder usually has 
no direct influence on the internal degrees of freedom, 
e.g. interaction among fermions. It is thus no wonder 
that the binding energy of the composite bosons exhibits 
no pronounced change, so does the chemical potential. 
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FIG. 2. The condensate fraction ric as function of (kpa)"^ for 
various disorder strengths. The shaded area and the vertical 
orange line at (kpa)'^ = represent the crossover region and 
the unitary limit, respectively. 



Condensate Fraction 

Though A and fi are the first quantities to study the 
crossover through the mean field theory but they are 
mostly of academic interest. From now on we shall focus 
on the quantities which are more prone to the experi- 
mental observation. Our first choice is the condensate 
fraction which is also one of our main results (depicted 
in Figj2]). In a clean Fermi gas, it is possible to work out 
the condensate fraction through mean field theory [44] 
which shows good agreement with the experiment |33j . 
Here we follow the similar mean field description, 



E 



r A(77) 



(13) 



where Uc is the condensate fraction. The only difference 
with the clean system calculation, in Eq.(13) is that, we 



used the disorder induced values for A and ii. We observe 
quite remarkable behavior for ric- As one expects the con- 
densate fraction decays similar to that of the clean limit 
when l/{kpa) — ?> —oo since A does not change with the 
variation of rj in the BCS regime. If one extends the self 
energy upto the second order in the condensate fraction 
calculation one can as well observe the effect of disorder 
in the BCS limit, where — n^g gets saturated at some 
finite value instead of exponential decay [55] , with de- 
noting the condensate fraction in the clean limit. In the 
EEC side, the disorder destroys part of the condensate 
and turns it into a normal fluid. The condensate fraction 
approaches roughly rj/y/kpa as obtained from the study 
of hard sphere Bose gas in random disorder [43] . 

The nonmonotonic behavior in the crossover region 
(grey area in Fig. ^ is the most intriguing point. In 



the study of quantized vortices, one sees that accumu- 
lation of a number of vortices becomes maximum in the 
crossover region [35]. Later on, in a theoretical study of 
a single vortex it was also observed that the circulation 
current is maximum in the crossover regime [35, . Sim- 
ilar feature was also reported in the Josephson current 
study |36j . This unique behavior was qualitatively at- 
tributed to the maximization of Landau critical velocity 
in this region. In general all these observations actually 
points to the robustness of the superfluidity at the uni- 
tarity. But normal mean field does not show any precise 
nonmonotonic behavior for the condensate fraction. Here 
with a weak disorder, we are able to generate a picture 
which affirms the belief that in presence of low amount 
of impurity the condensate fraction is less affected across 
the crossover, implying a comparatively high yielding of 
superfluidity in this region. 

One should also take a note of the position of the ex- 
trema. With the increase in disorder, we observe that the 
extrema slowly move towards 1/a — )• 0+ from +oo, but 
never cross it as we are exhausted with weak impurity 
limit. Further increase of disorder will break down the 
weakness condition. This result agrees qualitatively with 
Ref. [m ini] , however in those cases it has been suggested 
that the most robust region of superfluidity emerges when 
1/a — >■ 0_. But there is no good justification for that. 
Here it is pointed, this behavior is consistent with that 
of the critical velocity discussed in the following section. 



Sound Mode and Critical Velocity 

Another important and experimentally relevant |34j 
quantity is the lowest energy collective excitation of the 
condensate. The nature of sound velocity, in presence 
of disorder, in a Bose gas has been studied quite exten- 
sively, but it lacks a real consensus so that considerable 
amount of ambiguity still exists. Using a perturbative 
method, the sound velocity is enhanced in the presence of 
uncorrelated disorder and very weak interaction '161 1461 - 
I48j , whereas within a non-perturbative self consistent ap- 
proach and spatially correlated weak disorder case, de- 
pression in sound velocity has been reported [15H52) . In 
a more recent study it has been shown that there exists 
no generic behavior of sound in the presence of disor- 
der using a perturbative approach [53]. Hence in the 
crossover region, the behavior of sound is expected to be 
quite interesting. 

From technical point of view in order to obtain the 
sound velocity one needs to carry out an analytic continu- 
ation of the Matsubara frequency. Hence the fluctuation 
propagator matrix Ai is expanded to the second order 
for both momentum and frequency. The determinant of 
which leads to [54 , 



A^ii(q, i/)A^i2(q, -v) - A^L(q, ^) = 
^(A, li, k)q2 + B(A, ^i, k)iy2 + . . . ^ 0, 



(14) 



for v 



|q| where Vs represents the sound velocity. A 
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FIG. 3. The sound velocity Vs divided by the Fermi velocity 
vf as a function of (kpa)'^ for various disorder strengths. 
The shaded area depicts crossover region. 



and B are functions of A and fi and can be evaluated 
by summing them over k. In the BEC limit it has been 
shown that the sound velocity is — A^/(8m|/Lt|) [37] 
so that the suppression of the order parameter should 
directly result in that of the sound velocity. Near the 
unitarity, the sound velocity is directly connected to the 
chemical potential through = 2/i/(3m). In effect, the 
sound velocities for various disorders are merged near 
the unitarity as the chemical potential exhibits almost no 
change regardless of addit ion of the disorder. Moreover, 
in the BCS side, Vg = k"^ / {"iw?) ~ 0.57 (in dimen- 
sionless units) does not have any explicit dependence on 
the mean field parameters. Therefore, irrespective of the 
disorder the sound velocity saturates near 0.57. All of 
these arguments successfully explain the behavior of the 
sound velocity shown in Fig. [3] 

In an usual way, sound velocity is defined as Vg = 
where p is the mass density and kt is the isother- 
mal compressibility. Applying Gibbs-Duhem relation 
sound velocity can be written in usual notation as Vs = 

\ —'a^- An involved study of sound mode from the thcr- 

y m an ^ 

modynamical point of view reveals that the decrease in 
the collective excitations can be attributed to the in- 
crease in the compressibility. This is in accordance with 
Rcf. [49] . Intuitively this looks more feasible as the occur- 
rence of additional random potential should lead to addi- 
tional scattering and, resulting in the decrease of sound 
velocity [35]. 

Though the sound velocity itself is a quantity of huge 
interest, but it also serves an additional important infor- 
mation; According to Landau criterion, it determines the 
critical velocity of BEC. The critical velocity is defined as 
the minimum velocity of the BEC required to break su- 
perfiuidity by creating elementary excitations. To obtain 



FIG. 4. The sound velocity Vs (filled marks and solid curves) 
already shown in Fig. [s] and the critical velocity Vp (open 
marks and solid curves) obtained from pair breaking, normal- 
ized by the Fermi velocity vf, as a function of (/sfo)"^ for 
various disorder strengths. The critical velocity Uc's (solid 
curves) are determined from min{vs,Vp} . The unitary re- 
gion is depicted through the grey area. The inset focuses the 
turning point in the unitary region. 



a clear picture leading to the critical velocity, we calculate 
the minimum velocity related to the single particle excita- 
tions induced by pair breaking, which is dominant in the 
BCS side. The fermionic single particle excitation within 

mean field is written as {^J /S? + — /i) /m [37] , which 
is reduced to the familiar result of A/A; i? in the deep BCS 
region. In the case of the composite bosons, the Landau 
critical velocity is provided as the sound velocity. By 
choosing the minimum between Vs and Up, one can deter- 
mine the Landau critical velocity, which is represented 
by the solid curves in Fig. [4] 

Though the simple mean field analysis on Landau crit- 
ical velocity may not be amenable for a direct mapping 
between ric and Vc, but it may provide a qualitative sug- 
gestion for less depletion of condensate fraction in the 
crossover regime. In both cases the maxima exist near 
the unitary regime and move toward {kpo)^^ — Q as -q 
increases as shown in the inset of Fig. [4] 

IV. CONCLUSION 

In conclusion, we have studied several important physi- 
cal quantities such as gap parameter, condensate fraction 
and sound velocity etc. to address the issue of BCS to 
BEC crossover in a disordered environment. To this end 
we have included weak disorder via the gaussian fluctu- 
ation as prescribed earlier [26l [31] , and hence solved the 
coupled BCS mean field equations self consistently. This 
enables us to obtain the two basic mean field parame- 
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ters, A and fi. We thus show that the order parameter 
gets depleted leaving the chemical potential unchanged 
as we go from a BCS to a BEC regime. Afterwards, the 
condensate fraction has been calculated using the well 
known mean field description with the disorder affected 
A and As the disorder strength increases, we observe 
pronounced maximum developed in the crossover region 
justifying the expectation of robust superfluid in this re- 
gion. We have tried to connect this non monotonic nature 
qualitatively to the Landau critical velocity, which also 
shows a sharp maximum near unitarity. With increase of 
the impurity this peak slightly moves towards the unitary 
point. A similar feature is observed for the condensate 
fraction. In addition, the depression of the sound velocity 
is also been addressed, which might be related to the en- 
hanced scattering from the random scatterers employed 
in this model. 

To be precise, the the nature of A, Uc, and Vc is been 
reported here, when subjected to a weak random impu- 



rity. We hope our present study has shed some lights on 
physics of the BCS-BEC crossover in disordered systems. 
An interesting future perspective can be generalization of 
the theory for arbitrarily strong disorder and interaction. 
We also hope that these results would be observed and 
verified in experiment soon. 
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